On dispersive energy transport and relaxation in the hopping regime 



0. Bleibaum, H. Bottger 
Institut fur Theoretische Physik, Otto-von-Guericke-Universitdt 
Magdeburg, Germany 

V. V. Bryksin, A. N. Samukhin 
A. F. Physico-Technical Institute, St. Petersburg, Russia 
(February 1, 2008) 

A new method for investigating relaxation phenomena for charge carriers hopping between lo- 
calized tail states has been developed. It allows us to consider both charge and energy dispersive 
transport. The method is based on the idea of quasi-elasticity: the typical energy loss during a hop 
is much less than all other characteristic energies. We have investigated two models with different 
density of states energy dependencies with our method. In general, we have found that the mo- 
tion of a packet in energy space is affected by two competing tendencies. First, there is a packet 
broadening, i.e. the dispersive energy transport. Second, there is a narrowing of the packet, if the 
density of states is depleting with decreasing energy. It is the interplay of these two tendencies that 
determines the overall evolution. If the density of states is constant, only broadening exists. In this 
case a packet in energy space evolves into Gaussian one, moving with constant drift velocity and 
mean square deviation increasing linearly in time. If the density of states depletes exponentially 
with decreasing energy, the motion of the packet tremendously slows down with time. For large 
times the mean square deviation of the packet becomes constant, so that the motion of the packet 
is "soliton-like" . 
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I. INTRODUCTION 

In recent years much attention has been devoted to the 
study of relaxation processes of non-equillibrium charge 
carriers in strongly localized systems, where transport 
proceeds via phonon-assisted hopping, like photoexited 
charge carriers in band tails (see e.g. jjj- ||) and An- 
derson insulators (see e.g. 0- ||). In such systems par- 
ticular small relaxation speeds are observed. Often the 
smallness of the relaxation speed is attributed to inter- 
action effects. However, also in strongly localized non- 
interacting electron systems long lasting relaxation pro- 
cesses are known to be the rule, and not the exception 

The theoretical investigation of such relaxation pro- 
cesses is notoriously difficult, since the system is strongly 
disordered and always in a transient state. In most prob- 
lems one is interested in time scales which are large as 
compared to the time a single hop needs. For such time 
scales all quantities depend strongly on frequency, also 
for very low frequencies, so that the consideration of dis- 
persive transport is vital. On the other hand in most 
problems of interest both spatial and energetical disorder 
exists. Due to the latter fact the transitions are inelastic. 
The inelastic character of the transitions, the relaxation, 
leads to a flow of energy from the electron system to the 
phonon system. Due to disorder this transport of energy 
is also dispersive. Therefore, dispersive energy transport 
and relaxation are intimately connected. The investiga- 
tion of the relaxation process requires both the considera- 
tion of dispersive particle transport and the consideration 



of dispersive energy transport. 

The intricate physical situation manifests also in the 
equations, which have to be considered. Since the trans- 
port is inelastic already in the simplest approximation 
the investigation requires to find solutions to integral 
equations (see e.g. [jlO]]- [03). Due to the fact that the 
particle moves in an energy dependent density of states 
the kernel of these integral equations does not depend 
on the difference between the site energies only. The 
traditional method to handle the situation, the percola- 
tion theory, is not applicable here. The effective-medium 
methods by Movaghar and coworkers give, as pointed 
out by Movaghar and coworkers, wrong results for sys- 
tems at low temperature jf3| Q. Furthermore, beyond 
the Markovian approximation the derivation of these in- 
tegral equations represents itself an intricate problem. 
To our knowledge, beside the attempts by Movaghar and 
coworkers, numerical investigations (see e.g. H], [^)), and 
physical intuitive considerations (see e.g. [pi,), mainly 
Markovian equations have been used (see e.g. 

It is the aim of the present paper to fill this gap and 
to provide a formalism that can be used for studying 
relaxation phenoma of strongly localized, charge carri- 
ers far from equillibrium in the hopping regime taking 
into account both dispersive particle transport and dis- 
persive energy transport. To this end we focus, for the 
sake of definitness, on relaxation of photoxited, non- 
interacting, non-equillibrium charge carrires in band tails 
of, e.g. amorphous semiconductors, like amorphous Si:H. 
On the first glance, this problem seems to be rather spe- 
cial. If the number of charge carriers exited is small 
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Fermi-correlation is negligible, so that one has to cope 
with linear rate equations. The linearity of the transport 
equations is, of course, the basic ingredient in solving 
the problem. A closer look on the problems of inter- 
est reveals, however, that most of the problems can be 
formulated in this way. Clearly, how to achieve lineariza- 
tion depends on the experiment chosen, but, provided 
the number of charge carrires exited is small, always the 
smallness of the particle number can be invoked. In this 
case, the number of particles at a site is certainly not 
small as compared to its equillibrium value, but small as 
compared to unity. After linearization the structure is, 
in principle, quite similar as for the case considered here. 

Below we present the derivation of a framework for the 
consideration of relaxation phenomena due to phonon- 
assisted hopping at zero temperature. We focus on the 
limit of strong localization, where dispersive effects are 
expected to be most pronounced. Here the disorder man- 
ifests in a strong dependence of the transport coefficients 
on frequency for very low frequencies. Consequently, in 
this regime the diffusion propagator can not be calculated 
from Markovian transport equations, like those used e.g. 
in Ref . . To simplify our integral equations we use the 
concept of quasi-elasticity, the particle changes its energy 
only by a small amount by one hop. 

The method is applied to the study of dispersive en- 
ergy transport. We find that always two tendencies are 
present, if only the density of states not decreases with 
increasing energy. First, there is a widening of the packet 
due to the statistical spreading, the dispersive energy 
transport. Second, there is narrowing of the packet due 
to decrease of the density of states with decreasing en- 
ergy. The overall evolution is determined by the interplay 
of these two tendencies. We have studied both tendencies 
in two limiting cases, for a particle moving in a constant 
density of states and for a particle moving in an exponen- 
tial density of states. In the first case the impact of the 
density of states decrease is absent, so that only the sta- 
tistical spreading is present. Here the packet evolves into 
a Gaussian packet moving with constant drift velocity in 
energy space. The packet width increases with time as 
\ft. The other result is obtained for exponential density 
of states. Here both tendencies are present. We find, that 
the velocity of the packet is strongly slowing down with 
time. In that case the mean square deviation of energy 
eventually becomes time- independent. Consequently, the 
motion of the packet in energy space is "soliton-like" . 
The concrete results on the diffusion propagator, its time 
dependence and its moments for exponential densities of 
states are of relevance for photoluminiscense experiments 
on amorphous Si:H. 



II. BASIC EQUATIONS 

We consider photoexited, localized charge carriers in 
band tails at zero temperture. After exitation the charge 



carriers lower their energy by phonon-assisted hops be- 
tween localized states. Since T = 0, only hops from 
higher to lower energy occur. In this situation the charge 
carriers are in strong non-equilibrium. We assume, that 
the number of excited charge carriers is small, and their 
energies are far from the Fermi level, so that it is very 
unlikely that an electron jumps to a site already occu- 
pied. Consequently, we can neglect Fermi-correlation. In 
this case the electron transport can be described by the 
rate equation |l6| , |l7|| . 

~jr = ^2 (Pm'Wm'm ~ VmWmm') ■ (1) 
m' 

In calculating the transition probabilities we assume, 
that the electron-phonon coupling strength is weak, so 
that only one-phonon processes have to be taken into 
account. Then, at zero temperature, the transition prob- 
abilities are given by 

Wm' m = 9 (u) - E m < + £ m ) © (fm' ~ Em) W (|R mm '|) , 

(2) 

where 

W(|R mw /|) = i/e- 2a l R ™™'l. (3) 

Here a is the inverse of localized state radius, and v is 
the phonon frequency. The energy uj is the upper bound 
for the energy transfer by one hop. Note, that in the ma- 
terials of interest u> can be much smaller as the Debye- 
frequency, since not all phonons can interact with local- 
ized electrons equally well. Short wave-length phonons 
are ineffective since the electron-phonon coupling con- 
stant tends to zero for momenta q with q/2a ^> 1. There- 
fore, the effective upper phonon momentum is of the or- 
der 2a, and not of the order of the inverse lattice constant 
of the host material. Furthermore, in disordered systems 
the high-energetic phonons are localized, and need not 
contribute to transport by necessity. 

The first step function in front of the transition prob- 
abilities restricts the transitions to transitions between 
sites separated by at most u> in energy space. Thus, it 
decreases the energy relaxation speed. In impurity con- 
duction, and in nearly all papers on relaxation of charge 
carriers in band tails, this step function is usually re- 
placed by unity. In impurity conduction this is quite 
reasonable, since it is assumed that hops are restricted 
within narrow stripe near the Fermi level, which is small 
compared to Debye energy. In the band-tail problem, 
however, we can see no reason for neglecting it in ad- 
vance. 

To calculate the transport quantities of interest, we 
have to calculate the diffusion propagator. In order to 
render the analytical calculations feasible we introduce 
continuous coordinates. The change of representation is 
defined by: 

n (p) =^2PmS(p~ Pm) , (4) 
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where p = (R, e) and p m = (R m ,e m ). In this represen- 
tation the rate equation takes the form: 

^ = J dp'n(p')V(p',p), (5) 
where V is determined by the equations: 

V {p', p)= dp!T] (pi) w pi (p', p) , (6) 



«j w (p', p) = W (p', pi) [5 (px - p )-5tf- p)] , (7) 



where V is the volume of the system. The application of 
the averaging procedure to the structural factor rj serves, 
in particular, as a definition for the density of states, i.e. 

N{e) = { V {p)) . (15) 

Products of the structural factor rj are averaged accord- 
ing to: 

{V (pi) ■ ■ ■ V (Pn)) = M" (pi) 6 (pi - p 2 ) • • • 6 (p n _i - p n ) . 

(16) 



v (p) = Y / Hp-p m )- (8) 

m 

The Laplace-transformed equation is given by: 

sn (p) - n (p) - J dp'n (p') V (p', p) , (9) 

where n Q (p) = n(p, i = 0) is the initial condition. Wc 
will assume that p m (t = 0) is a function po(R m , s m ), so 
that n a (p) = pq(p)t](p). 

Equation (^) can be solved using the Green's function 
method. The solution is given by: 



P = 



+ O 



o 



n(p) = J dp'n (p')<P(p',p). (10) 

The Green's function satisfies the equation: 

S $(p',p)- J dp 1 V(pi,p 1 )<t>( Pl ,p) = S(p>-p). (11) 

Note, that, due to probability conservation, the Green's 
function $ and the probability Wp satisfy the relations: 

J dp$(p',p) = i, (12) 



dpw p (p',p) =0. (13) 



III. CONFIGURATION AVERAGE 



FIG. 1. Diagrams contributing to P. The leading empty 
point represents the single structur factor rj. Everey full dot 
is associated with a potential V{pi,pj). Consecutive points 
are connetcted by <5(p; — Pj)/ S , represented by a solid line. 
Integration is performed over intermediate arguments. 
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FIG. 2. Diagrams contributing to the confguration average 
of P. Averaging is symbolized by dashed lines. Correlated 
averages are symbolized by joint lines. 



In order to calculate the configuration average, we as- 
sume, that the sites are distributed homogeneously in 
space. The distribution of site energies is supposed to be 
given by some distribution function p({ei}). Accordingly, 
the average of any quantity, depending on the energy and 
positions of sites, is given by: 



(A) 



n„ 



■de m p({e m }) A({R m ,e m }) , (14) 



n = 4 



-* + 



FIG 3. Diagrams contributing to II. Here, in contrast to 
Fig.jj], every full line is associated with a function F(pi, pj). 
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FIG 4. Diagrams contributing to S. Here every full line is 
associated with a function F(pi,pj). 
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Using these definitions, the configuration average of 
Eq. ( |Io|) can be calculated diagrammatically JlJJl!^] . The 
diagrammatic method leads to the following set of equa- 
tions for the calculation of the configuration average 
(n(p)) of the electron density n(p) pjfl , (see Figs |l|-^ 
for illustration): 

( n (p))= / dpidp2Po(pi)S(pi,p2)F(p 2l p) , (17) 



sF (p', p) = S(p'-p) + l d Pl n (//, Pl ) F ( Pl ,p) , (18) 



Tl(p',p) = J d Pl M( Pl )IL pi (p',p) 



(19) 



n pi (p',P) = w pi (p',p) 

dp 2 dp 3 w pi (p', p 2 ) F (p 2 , P 3) n pi (p 3 , p) , (20) 



SV,p)=^( £ ') 



5{p -p)+ j d Pl F (p, Pl )U pl ( Pl ,p) 

(21) 



Here F(p',p) = ($(p',p)). A detailed derivation of the 
set of integral equations is given in Appendix ^[ 



IV. EFFECTIVE MEDIUM 



Given the system of integral equations (|lj)-(21), the 
main problem is to find an approximate self-consistent 
solution to it. The situation is quite similar to that of 
the calculation of the equilibrium conductivity in a disor- 
dered system [[l8],[l9]. There an approximate solution of 
this system could be found by introducing a proper de- 
composition of the function F, the diffusion propagator, 
into short- and long-wavelength limit, according to: 



F{p>,p) = f{ s )C{e)5{p l -p) + F(pi,p) , 



(22) 



where f(s) was a frequency dependent parameter that 
could be related to the critical hopping length R c via the 
equation: 



fv = exp {—2aR c 



(23) 



and C{e) was an energy dependent function determined 
by the principle of detailed balance. This decompostion 
originated from the notion, that the integrals in the in- 
tegral equations are governed by products of transition 
probabilities and diffusion propagators, and, for strongly 
localized electrons, the latter quantities are short ranged 
functions as compared to the transition probabilities. Us- 
ing the decomposition (p2|), the effective medium approx- 
imation reduces to the repacement of F by fC(e)S(p' — p) 



in the calculation of the effective transition probability n 
(Eq.©) and of the irreducible block S (Eq.©). 

Here in the band-tail problem we use the same philoso- 
phy. We first decompose the diffusion propagator in two 
parts, according to: 



F(p',p) = f(e, S )5(p'-p)+F(p',p) 



(24) 



Then, to investigate the consequences of this decomposi- 
tion, we insert Eq.(|24|) into Eq.(po|). Performing a partial 
summation wc obtain: 



n p {p\p) = w~p (p',p) 
+ / dp 1 dp 2 Wp(p,pi)F(p 1 ,p 2 )Tlp(p2,p), (25) 



where the renormalized transition probabilities Wp(p' ,p) 
are given by Eq.(Q), with W replaced by: 



W(p',p) 



l + f{e',s)W{(/, P ) 



(26) 



Note that the renormalized transition probability de- 
pends now on s via the function /(e, s). 

At this point the introduction of the function /(e, s), 
the effective medium, still seems to be rather arbitrary. 
However, if we now choose the effective medium in such 
a way, that the integrals over F vanish, the renormal- 
ized transition probabilities yield the exact solution of 
the diffusion problem. In that case the function S turns 
into: 



S(p',p)=Af(e')6( P '- P ) 



so that Eq.(^) can be cast into the form: 



(n(p;s))= / dp (n (p)) F (p,p) 



(27) 



(28) 



Therefore the function F( P ' ', p) can be identified with the 
diffusion propagator. 

Note, that although the general philosophy is quite the 
same as in equilibrium, the situation is much more intri- 
cate here, — the principle of detailed balance is absent, 
since we are dealing with the situation at zero tempera- 
ture, and so it can't determine the energy dependence of 
the effective medium. Moreover, we expect the critical 
hopping length to depend somehow on the position of the 
electron in the tail. Thus, in contrast to equilibrium, / 
is not only a frequency dependent parameter but also an 
energy dependent function, which has to be determined 
self-consistently. 
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V. DIFFUSION PROPAGATOR IN EFFECTIVE 
MEDIUM APPROXIMATION 

For the moment we put aside the question of the de- 
termination of the effective medium to elaborate further 
on the consequences of the renormalization of the transi- 
tion probabilities. To this end we focus on the diffusion 
propagator. 



The equation for the diffusion propagator is given by 
Eq. ( |l4[ ) . In effective medium approximation, when calcu- 
lating the irreducible part II, only the lowest order con- 
tribution to II with respect to F is taken into account, 
so that lip = Wp. In that approximation the equation 
for the diffusion propagator in momentum representation 
reads: 



sF (qje 7 , e) = S{e' -e) + J efci {f (q|e', e x ) W (q|ei, e; s) N ' (e) - F (q|e', e) W (0|e, £ X ; s) AT( El )} 



(29) 



Of course, as we don't know how the effective medium looks like so far, and moreover, as the equation is a complicated 
integral equation, we can't find a solution. The fact, that the effective medium is a function of energy, makes the 
problem much more complicated. Progress can only by achieved if we can find arguments to simplify the equation 
considerably. To this end we focus on the renormalized transition probability. According to the Eqs. (|^) and (p6[), it 
is given by: 



W (R\e', e;s) = Q (V - e) 9 (a - e' + e) 



W (R) 



l + f(e',s)W(R) 



= 6 (e' - e) Q (w - e' + e) W (R\e'; s) 



(30) 



Owing to the step functions in front of the transition probability, the energy integrations are restricted to intervals of 
length bj. Taking this fact into account the integral equation for the calculation of the diffusion propagator takes the 
form: 

sF(q|e/,e) = 6(e'-e) + J tfei {f (q|e', e + e x ) W (q|ei + e; s) Af (e) - F (q|e', e) W (0|e; s) Af (-e x + e)} . (31) 



Further we assume that the diffusion propagator, the transition probability W(R,e; s), and the density of states 
jV(e) are slowly varying functions on intervals of length w, so that the integrand can be expanded with respect to e±. 
Doing so we obtain: 



sF(q\e/,e) =8{e' -e)+wM{e) \w (q|e; s) - W (0|e; s)] F (q|e', e) + \^ 2 ^- £ 



F( q \e , ,e)W(q\e;s)N(e) . (32) 



We terminate the expansion after the first derivative 
with respect to energy. This term describes the biased 
motion of the particle from sites of higher to sites of 
lower energy. For finite temperatures one would have 
to replace F(q|e',£) in the last term with F(q|e',e) + 
kTdF (q\e' , e) jde which would describe the energy dif- 
fusion current. At zero temperature there is no input of 
energy from the phonon-system into the particle. There- 
fore, we expect the additional term to be rather small, 
and can be neglected. 

To investigate the transport properties, we restrict our 
attention with the long- wavelength limit. In that case 
the elastic contribution yields the diffusion coefficient for 
particle diffusion. The term, containing the derivative 
with respect to energy is finite for q — > 0. If we are inter- 
ested only in the long- wavelength limit, we can put here 
q = 0, since the remaining terms are of higher order with 
respect to 'g' and , u>d/de\ Since a nonzero momentum in 
this term couples particle diffusion to energy transport, 



this approximation corresponds to a decoupling of these 
two processes. Then, in the long-wavelength limit, we 
obtain: 



sF (q|e', e) = 5 (e' — e) — D (e, s) q 2 F (q|e', e) 
d 



where: 



and: 



— [F(q\e',e)v(£,s) 



D(e,s)= -~wV 2 q W(q\s;s) 



q=0 



Eq.(tS3J) is easily solved. It's solution is: 



(33) 



(34) 



(35) 
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F(q| £ / ,e) = 



9 (e' - e) 
v (e; s) 



■ exp 



de\- 



-D{ei\s)q 2 
v{ei\s) 



(36) 



Explicit expressions for the transport coefficients can be 
found in Appendix 

Now, at this stage, the validity of our quasi-elastic ap- 
proximation requires, that the second derivative terms 
are small as compared to terms with first derivatives with 
respect to energy. This requirement imposes the follow- 
ing restrictions on the transport coefficients v and D, the 
frequency s and the momentum q: 



1 dv 

v de 



«1, 



«1, 



2 D , 

q z uj— < f . 

v 



(37) 



(38) 



(39) 



The applicability of th e q uasi-elastic approximation was 
also discussed in Ref. |p0|| . There it was concluded that 
this approximation should be inapplicable. To substan- 
tiate this statement numerical calculations were invoked. 
However, in interpreting these data it has to be taken into 
account, that they have been obtained using a model, 
neglecting a weighting of the transition probabilities ac- 
cording to the number of phonons emitted, so that hops 
between nearly isoenergetic sites were treated as likely 
as hops from the very top to the very bottom of the 
tail. Consequently the discussion in Ref. [^o| applies only 
to system with sufficient strong electron-phonon interac- 
tion, but not to systems with weak electron-phonon in- 
teraction. For extraordinary deep hops to contribute to 
the diffusion propagator they should be characteristic for 
an ensemble of electrons. This is, however, not expected 
and in the experimental data, e.g. on amorphous Si:H 
pjj , not observed. For these reasons the discussion in 
Ref. pol does not apply to our model. 



VI. PHYSICAL QUANTITIES AND 
INTERPRETATION 

Before establishing a self-consistency equation, we first 
elaborate further on the physical content of our diffu- 
sion equation. To this end let us first have a closer look 
on the coefficients D(s,s) and i>(e, s). Imagine, that we 
have a particle initially located at (Ro,£o)- Then the 
initial condition is (n (R, e)) — 5(R — R )<5(e — £o)- Ac- 
cording to our approximation, the motion of the charge 
carrier is composed of two contributions: particle dif- 
fusion between isoenergetic sites, and relaxation in en- 
ergy space. Characteristics of these two processes are 



the mean square displacement and the energy relaxation 
speed. It turns out, that both can be calculated from the 
function: 



u t \ © (e - e) d 
s ae 

6 (s - e) 



de\ 



v(ei,s) 



v (e, s) 



■ exp 



efei 



v (ei, s) 



(40) 



which is just F(q = 0|eo,e)- 

We define the mean energy of the particle by: 

E(s)= [ dp'e' (n(p',s)) . 



(41) 



Using our diffusion propagator, this equation can be 
rewritten in the form: 



E{s) = J de'P L {e ,e';s)e' , 



(42) 



and from the diffusion equation we obtain after integra- 
tion by parts: 



sE (s) - e = - J deP L (e , e; s) v (e, s) 



(43) 



Therefore, in general, the time dependence of the mean 
energy is given by complicated integrals. These integrals 
simplify considerably in two limiting cases: in the ab- 
sence of dispersive energy transport, and for energy in- 
dependent v(e,s). In the absence of dispersive energy 
transport, i.e. in the Markovian limit in which the trans- 
port coefficients are independent of s, the integrals can 
readly be calculated in time representation. In the latter 
situation we simply obtain: 



P L (e , e;t) = Q (e - e) 8 (e m (t) - e) , 
where e m {t) is defined by: 



* = 



de\ 
v(et) 



(44) 



(45) 



Therefore we obtain: 

dE (t) _ de m (t) 



dt 



dt 



= -v (e m (i)) 



(46) 



Consequently, v is the velocity of energy relaxation, and 
e m (t) is the instantenous position of the particle in energy 
space. In general, however, in disordered systems the co- 
efficient v(e, s) depends on s, so that energy transport is 
dispersive. In that case the integral in Eq.(^3|) can only 
be calculated easily if v is independent of energy. Then, 
owing to probability conservation, 
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dt / V ' s 



is obtained. 

The mean square displacement, defined by: 



(47) 



R z {t)= dp'R' 2 (n(p')) , 



(48) 



can be obtained from similar arguments. In Fourier rep- 
resentation this equation can be written in the form: 

^) (s) = sA q \ q=Q J de'F (q\e , e>) . (49) 

The derivative of the diffusion propagator can be calcu- 
lated using the diffusion equation. Doing so, we obtain: 

(If) (S) = 2d J d£Ph (£ °' 65 S) D (£ ' S) • (50) 

Again Eq.(pO[) simplifies only in two special cases, in the 
Markovian limit and for energy independent diffusion co- 
efficients. Whereas in the first case 



d ^(t) = 2dD{e m (*)) 



is obtained, we have 
'dR 2 



at I s 



(51) 



(52) 



in the latter situation. Below we shall see, that energy 
independent transport coefficients are obtained for con- 
stant density of states only. 



VII. THE SELF-CONSISTENCY EQUATION 

So far we have only investigated the consequences of 
the renormalization. In order to complete the approx- 
imation scheme, we still have to calculate the effective 
medium itself. Of course, we can not calculate the effec- 
tive medium exactly, since this would amount to find an 
exact solution to the diffusion problem. Rather we shall 
try to calculate f(e, s) self-consistently. 

The transport coefficients are properties of the Green's 
function F, the diffusion propagator. Thus, in order to 
find an equation for f(e,s) we should relate the trans- 
port coefficients to II, the irreducible part of the diffu- 
sion propagator. The diffusion coefficient comprises only 
elastic contributions, however, in the relaxation problem 
a proper description of inelastic processes is vital. More- 
over, as shown in the previuos section, the character- 
istics for particle and energy transport can already be 
calculated from the function Pj,, the q — > limit of the 
diffusion propagator. Therefore, in establishing a self- 
consistency equation we should focus on those character- 
istics that are important in that limit, that is the energy 
relaxation speed v. 



The equation for the diffusion propagator for q = is 
given by: 

sP L (e',e; s) = 5 (e' - e) + P L (e',e, s) J deiU (si, e; s) 
dP L (e',e; s 



ds 



dei (ei — e) II {s\,e] s) 



(53) 



If we compare Eq.(^) with Eq.(^) we deduce that 

v(e,s) = J de% (si — £)n(ei,£;s) . (54) 

Now we decompose v, as defined in Eq.(j54|), into two 
parts, one part that contains only the effective medium 
approximation of II: 

v(s,s) = J de 1 (e 1 - e) U\ EMA (e u e; s) , (55) 

and a part 6v(e, s), that contains the deviations F. Self- 
consistency requires: 



5v (e, s) = . 



(56) 



v(e), as defined by Eq.(p5|), is in accordance with the 
definition (^5|), taking into account Eq.(B3) and the in- 
equality Ljf(e)/f(e) < 1 (Eq. @), owing to which 
contributions proportional to this parameter are negligi- 
ble. 

We now focus on the self-consistency equation (|5q). 
While Eq.(|5|) contains only the effective medium contri- 
bution to the diffusion propagator, Sv is a functional of 
the effective medium / and the deviation F. By construc- 
tion, it is at least linear in F. If this equation could be 
solved exactly, an exact solution to the diffusion problem 
could be found within quasi-elastic accuracy. In practise 
this is not possible, and, therefore, we depend on fur- 
ther approximations. To simplify this equation, we take 
into account only the lowest order contributions to this 
equation with respect to F, i.e. we linearize dv with re- 
spect to F and require that the first order contribution 
vanishes. This approach is quite close to the usual CPA- 
philosophy, in which vanishing of the t-matrix is required 
in its lowest order approximation. Using this procedure 
we obtain the following self-consistency equation: 



1 



uj 2 Af(e)f(e,s)W(0\e-s) = auj 



sb 



W{0\e;s)Af(e) 



(57) 



where a and b are simply numbers. A detailed derivation 
of this equation is given in Appendix 

In deriving the self-consistency equations we have im- 
posed further restrictions on the effective transition prob- 
abilities, which determine the range of its applicability. 
These restrictions can be formulated most convinicntly 
using the dimensionless critical hopping length p c (e,s), 
related to f(e, s) by (see Appendix H): 
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/(e,s)^ = expp c (e, s) 



(58) 



In terms of p c , the inequalities ( |37| ) and (|38|), used in the 
derivation, read (prime is derivative with respect to e): 



\u>p' c (e,s)| « 1, 



(59) 



explicit equations for the transport coefficients can be 
obtained. They are given by: 



D (e, s) D (e, s) s 
In ■ — 



D(e,0) D(e,0) 0(e)' 



(65) 



\p c (s,0)-p c (e,s)\ «p c (e,0) . (60) 
In addition, when calculating the integrals, 

p c (e,s)>l, (61) 

was used. 

A closed solution to the self-consistency equation can 
only be found in the limit s = 0. There we obtain: 



p c (e,s= 0) = 



2a 



l/d 



2da 



S(d) 



l/d 



(62) 



where S(d) is the solid angle in d dimensions. 

For s satisfying (pCf) , Eq. (R7H can be cast into the form: 



[p c (e, 0) - p c (e, s)] exp [p c (e, 0) - p c (e, s)] = 



0(e) 



where: 



0(e) 



6w/9 c (e, 0) 



v(e,0) . 



(63) 



(64) 



According to Eq. (|63|) , the critical hopping length de- 
creases with increasing frequency. Note, that the struc- 
ture of the equation (|6^) for the calculation of the disper- 
sion of the critical hopping length, obtained here, agrees 
completely with that obtained for the critical h opp ing 
length in calculating the equilibrium conductivity |]l6[|l9| . 

The dispersion of the transport coeffcicients is deter- 
mined completely by the dispersion of the critical hop- 
ping length. For small s the frequency dependence pre- 
exponetial factors can be ignored, so that from Eq.(p3) 



v (e, s) v (e, s) s 
v{e,0) n v(s,0) = fi(e) 



(66) 



The formal solution of these equations is given by the 
Lambert's W-function W(z), defined by the equation 
z = W(z) exp W(z). Using the Lambert's function we 
can write: 

D(e,s) = J D(e,0)expW(s/^(e)) , (67) 



v(e,s) = v(s,0)expW(s/n(e)) . (68) 



VIII. CONSTANT DENSITY OF STATES 



A constant density of states, although of not much 
physical relevance, gives us the unique opportiunity to 
study pure dispersive energy transport. Here both O(e) 
and wo = v(s,0) are independent of energy. Conse- 
quently, Eq.([40|) can readly be integrated. The integra- 
tion yields: 



Pl (eo -e,s) = -— — exp 

v[s) 



s (e - e) 



v(s) 



(69) 



The time-dependence of this function can be obtained 
by the inverse Laplace-transformation. Using Eq.(|6^) to 
change the integration variable from s to y = v/vq, the 
inverse Laplace-transform of Eq.(69) may be written as: 



P L (e -e,t) = — 0(e o -e) / 
«o Jc 



dy In y + 1 



2iri 



1J 



■ exp 



Qty In y - — (e - e) In y 



(70) 



with the properly chosen integration contour C. At large enough t, that is for fit 3> 1, this expression, using the 
saddle-point method, gives simply the Gaussian packet: 



Pl {so - e, t) w 6 (e - e) 



n 



i 



v 2^-K^lt 



exp 



(e -s~ v t) 



(71) 



Now let us consider the time dependence of the energy dE(t) _ j ds v (s) ^ st 

relaxation. According to Eq.(^) the velocity of energy dt J c 2ni s 

relaxation is: 

Again, the time dependence of the energy relaxation 
speed can be calculated using asymptotics. Details of 
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the calculations are presented in Appendix [S], where it is 
shown that: 



dE (t) 
dt 



- Vq 



■ exp (— Clt/e) , as ilt ^> 1 



(nty 

^[\n(e/(nt))}-' 2 , asflK< 1 



(73) 



Note, that the problem of the energy relaxation in the 
case of constant density of states (v(e, s) is independent 
of e) is completely equivalent to one of the nonmarkovian 
charge transport in strong electric field E, when the diffu- 
sion (described by the second coordinate derivate) is to- 
tally neglected. One has only to replace v (s) — > u (s) E, 
u (s) being the mobility jl5],[l6) . In the context of disper- 
sive particle transport, the regimes fit <C 1 and fit ^> 1 
are the regimes of anomalous and normal "diffusion" , re- 
spectively. The main difference between both is in that 
while the sites are ususally distributed homogeneously 
in space, the density of states is usually an increasing 
function of energy. 



IX. EXPONENTIAL DENSITY OF STATES 

A. The saddle-point approximation and its break 
down for large times 

The calculation of the time dependence of the energy 
distribution function for an arbitrary density of states on 
the basis of the Eqs. (§|) and © turns out to 

be a quite intricate problem. A tool that can be utilized 
in tackeling the problem is the saddle point approxima- 
tion. How to apply the saddle point approximation for 
the calculation of the quantities of interest for an arbi- 
trary density of states is shown in Appendix || Below 
we focus on the exponential density of states, which is 
relevant, e.g., for amorphous Si:H. 

We assume that the density of states is given by: 



Af(e) =A/- exp(3-l) . 



(74) 



To simplify the notations we use the abbreviations: 



CO 



buj 
2da 



(75) 



p(e) = p c (s,0) = Aexp(-^) , 



where: 



A 



2a 



oAf(Q)] 



l/d 



2da 



1 1 / d 



S{d) 



(79) 



(80) 



Using these equations it follows from the formulas de- 
rived in Appendix |e] that the time dependence of the 
mean squared deviation and the time dependence of the 
mean energy are given by 



ct 2 (e m , eo) = c 2 (i, e ) ~ wA 



1 - 



(*) 



where 



-A In 



to (eo) 



1 , . 



(t + toY 



to) 



vp (e ) v (e Q ) 



(81) 



(82) 



(83) 



Note, that, according to the saddle-point approxima- 
tion, the distribution is Gaussian. Furthermore, the dis- 
persion is constant for t 3> to- Consequently, for time 
scales in line with the applicability of the saddle-point ap- 
proximation, the motion of the energy packet is "soliton- 
like" , that is the packet moves without distortion. 

The applicability condition for the saddle-point ( E13| ) 
requires: 



P 2 {e m ) < 1. 



(84) 



Because of p(e m ) grows with t, the saddle-point method, 
and, consequently, all the results of this chapter, be- 
comes invalid, at least at sufficiently large t. The phys- 
ical meaning of the condition J84|) is the following: at 
t 3> to the width of the electron's energies distribution 
is Se — \AlA. Because of p ~ exp (— e/A), the variation 
of p is: Sp/p = y/ui/A, and, because v ~ exp(— p), we 
have: Sv/v = 5p = y/ujp 2 /A. Thus, the condition (|84|) 
is just the one of the variation of electron's "velocity" 
across the distribution is smaller, then the velocity itself. 



2da 2 



(76) 



Then, the Eqs.(g3j), (|B2|) and (|6J) can be cast into the 
form: 



' (e, 0) ee v (e) = ujv exp [~p (e)} , 



n( e ) 



exp [-p (e)] 



(77) 



(78) 



B. Form of the distribution for large times 

Even if the initial condition are such, that the con- 
dition ( jS4| ) is fulfilled, and a gaussian distribution is 
formed, at some moment of time p m = p(s m ) becomes 
of the order of \J 'A/a;, which definitely should results 
in some deviation of the distribution function from its 
symmetric, gaussian form. For the exponential density 
of states, distribution function in the Laplace represen- 
tation, Eq.(Eoh, using Eqs . (|68|,[77|-|79|) , may be written as: 
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P L (£o, e; s) = — p— exp 



fix 



W(sxe x ) 



Pa 



W jspeP ) 
uips 



■ exp 



dx 



W (sxe x ) 



Po 



(85) 



where it was set v = 1 and A = 1 by appropriate choice of time and energy units. Under the same conditions, which 
were used in the derivation of the above formula, the following approximation is valid (see Appendix [f]) : 



W<Hre x ) 1 
dx 1 p > ^G(s,p)~G (s, po) ; G(s, p) = -W (spe?) 
ar p z 



1 + -W (spe p ) 



(86) 



When calculating distribution function at sufficiently large times, the characteristic values of s, giving the main 
contribution into inverse Laplace integral become so small, that one can suppose |s| po exp(po) 1, and: 
G(s,po) w spQ 1 exp(po)- Under this condition the distribution function approach initial conditions independent 
shape: P(e , e; t) -> 0(e, t + to(eo)), to(e ) = ifiPof expp , 



+ 200 



4>(e,t) 



ds 



2m 



-0 L (e,s)e s 



-\-ioo 

f ds d f 1 ,,,, fl , 

/ 7T—PT r exp\st-—W(spe p ) 
J Ittis op i ojp z 



1 + -W(spe») 



(87) 



Let us consider, first, momenta of the distribution at the large times. In Laplace-representation the momenta are 
defined by the equation: 



oo oo 

Xul(s) = J dte- st (e n ) (t) = J ^(-lnp) n P L (e ,e,s)^e st ° j ^ 



(-]xip) n 4> L (e,s) 
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We shall omit further on irrelevant time shift multi- ~ _j_ y/fi) \ n 

pie e Bt °. It turns out, that for Cop ^> 1, the integrand of Pm ~ y J n 

Eq.(^||) has a maximum, located at: 



However, if one would try to use the saddle-point method 
in evaluation of integrals like (|8q), the results would be 
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incorrect at \[Qp m 3> 1, since the saddle point criterion 
is not fulfilled. However, at cup 2 ^s> 1, momenta can be 
calculated as an expansion on the powers of small para- 
menter l/(-\/15p). Details of the calculations are pre- 
sented in Appendix [f| 

As a result, we have for the first momentum, or the 
mean energy: 



Xi (t) = (e) (t) 



1 . (~ 
A 



A In — In (bVtit) 



1 



8w m 2 ( i 



(89) 



where b = \f2e~ 1+ ~</ 2 , and 7 is Eulers's constant. For the 
distribution's dispersion cr 2 (t) = ^{t) = X2(t) — XiCO we 
obtain: 



^(.)- MW -( J -|) + (l + laJ)^ + ..., 

(90) 

and for the third central momentum (J-a(t) — 
((e-<e) (t)) 3 ) = X3 (f)-3 X i (t) X2 (t) + 2%? (t) we have: 

M3 W = -( 7 r-3)y|<I ) 3 / 2 + 3[ 7 r(l-ln2)-l] i -^=^, 

(91) 



Figures |[ ^, and [?] show mean energy, xii mean square 
deviation of energy er, and dimensionless coefficient of 
asymmetry: 



Asymmetry = 



M3 (t) 
a* (t)' 



(92) 



respectively, plotted versus In Inf. There values, ob- 
tained numerically, using the distribution function (|87|), 
are compared with their corresponding asymptotic forms, 
following from Eqs. 



Moving along the same line, the whole distribution 
function may be reconstructed at large times, if to cal- 
culate all its momenta. However, it may be obtained in 
more simple way. Namely, when up 2 3> 1, one may re- 
place s within arguments of W- functions in Eq. ( |S7| ) with 
some so — c/t, c ~ 1 (a kind of argumentation may be 
found in Appendix |G]). Then the integral can easily be 
calculated to yield: 



(e,t) 
d 



p— exp . _ 
up 1 ujp z 



1 t< 



pe 



1 

2 W 



: pe' 



(93) 



a) 




Energy 



b) 




Energy 





Energy j ^ H ' _~ ' ' _ 4 ' ' _ 3 V ^ Energy 

FIG. 8. Plots of the energy distribution function at Q/A = 0.1 and: (a) t = 10 3 , (b) t = 10 4 , (c) t = 10 6 , and (d) t = 10 9 . 
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p 



p 




p 



p 




One can easily show (see Appendix |g|), that asymp- 
totic expressions for momenta ( p9[|9"l| ) may be entirely 
reproduced from Eq.(|93|), if to set c = exp(— 7) (7 is the 
Euler's constant). Figures [8[]l0| show some plots of the 
distribution function (87), compared with its asymptotic 
form (^3|). This latter can be further symplified, if to 
introduce pt: 



- = Pt e pt ; p t = W ( - 
c \c 



(94) 



As we shall see, the body of the distribution at a given 
time corresponds to values of p > p t , p — pt <C pt- Then, 
one can use approximation: W {(c/t) pe p ) ~ p — pt^> 1. 
Then, we have in the exponent (1 — p t j p) /2uj, which 
can be replaced by (e — e t ) /2u), where e t ~ lnlni. As a 
result, at large time and small Hi the distribution function 



turns into: 



(e-e t f 



2ui 

1 , e > e t 



< e t 



Et = In In t . 



(95) 



C. Initial stage of the evolution 



At t = it was supposed P(e,eo;t = 0) = <5(e — e ). 
To understand, what happens at the intitial stage of the 
distribution's evolution, let us consider the behaviour of 
P at e « £ Q . One can write from Eqs.(||J6|,0]79|): 
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P (e ,t) = P(e a ,e ;t) = -J- [ -J-W (~p eP°) 

ujpo J 2ms \v / 



1 



dz 



- (1 + z) exp (rze z ) 
1 f dz 

— / ^— exp (Tze - z) ; 
j t J c 2m 

po = A exp 



ujpo Jq 2ni 
1 

uipo 



m 



-po 



(96) 



where the new integration variable z = W {{s/v) poe Po ) 
was introduced instead of s, and the integration by parts 
was performed. This integral may be evaluated using 
the saddle-point method. The saddle-point equation, 
f'(z) = t(z + l)e z — 1 = 0, gives us saddle point value 
z c = W(e/r)-l, that is z c « ln(l/r) - In ln(l/r) at r < 
1, and 2 C « e/r — 1 as r > 1. To ensure the correctness 
of the saddle point approximation, one should require the 
parameter|/'"(z c )| 2 /l/"(^)| 3 = |(3 + z c ) 2 /(2 + z c ) 3 \ to 
be small. While this is true at t > 1, this parameter 
appears to be w 1 at t < 1. Therefore, the result for 
t <C 1 is correct up to the multiple of the order of 1 only. 
We have: 



(f) 3/2 exp(-f) as r> 1, 



Po(p,t) 



/2TTOp r V el ' f07") 

—t=^ — (in — — In In — ) as t 1 • 



In a similar fasion one can calculate the energy deriva- 
tive at e = £q: 



Px(po,t) 



1 



-t-oo 



LuApo 
1 

(vpof 



dP 
~de 

ds 
2nis 

dz r 
2Tri 



-W 



/AW W 



P{po,P\ t) 



z(l + z)- 



1 + W 



-Po 



exp(rze z ) , (98) 



where the term (uj/A)pqz within the square brackets was 
neglected due to assumption (tD/A) p -C 1. At t ^> 1 this 
integral may be readily evaluated in the saddle point ap- 
proximation, which yields: 



Pi(p,t) 



1 



2ttt 



exp 



(- 1 
V e 



(99) 



When r <C 1, one have to perform integration by parts 
first, which gives: 



Px{p,t) 

dz 
q 2ixi 



1 



UJ 2 p 2 T 



z-1- 



LUp 



ujp 



exp (rze z — z) 



A(l + z) A(l + z) 
Then one have in the saddle point approximation: 



2nu 2 p 2 



LUp~ 



In- 

T 



In In — . 

r 



(100) 



At sufficiently small t the above expression is positive, 
which means that the distribution is "sticked" near s = 
So, that is, it monotoneously decreases as eo — £ > 
increases. At some r = To, Pi changes its sign. If 
(Q/A)pq ^> 1, this corresponds to z c0 s» In — w 

^/w/Apo, otherwise tq ~ 1. This correspons just to the 
moment of time, when this distribution separates from 
the intial point, — its maximum is at e < eq. For r <C tq 
the width of the distribution on the energy scale may be 
estimated as: 



Ae(t) 



up 



Po{p,t) _ ujp 

Pi Go,*) ~ z c ~ ln(l/r) 



(101) 



This formula is valid if ln(l/r) <C p, up to the times cor- 
responding either r ~ 1 if (Co/A)p 2 t -^i 1, or ln(l/r) ss 
y/ui/Apo otherwise. In the former case the distribution 
width becomes of the order of uip just before it separates 
from the initial point. After this, the distribution, as 
it was shown in previous subsection, widens till y/u>A, 
and the gaussian packet moves downward, untill its cen- 
ter reaches the value, corresponding to p ~ y/ A/a). In 
the latter case the packet's width is ~ V&A Hip at 
the very moment of its "unsticking". In both cases, to 
consider subsequent evolution, one has to use some other 
approximation instead of saddle point one. 



X. CONCLUSIONS 

In the paper we have presented an effective method 
which pcrmitts the investigation of relaxation phenom- 
ena of localized charge carriers far from equilibrium due 
to phonon-assisted hopping at zero temperature. From 
the point of view of the formalism the main equations are 
the equations for the calculation of the diffusion propa- 
gator JP), @ and ©, and the equations (g), (||) 
and (|66|)7 which determine the dispersion of the trans- 
port coefficients. These equations show, that both par- 
ticle transport and energy transport are dispersive. The 
equations ( |65| ) and ( |66[ ) lead to a strong dependence of 
the diffusion constant D(e; s) and the velocity of energy 
relaxation w(e; s) on frequency s already for low frequen- 
cies. To our knowledge, these equations have not been 
dervied in the literature so far, for systems far from equi- 
librium. In fact, in the literature mainly frequency in- 
dependent transport coefficients can be found (see e.g. 
pfH and references therein). The strong dependence of 
the transport coefficients on s results in a non-Markovian 
equation for the calculation of the diffusion propagator 
(p3|). This discriminates our equations from the Marko- 
vian integral equations used in the Refs. Jjc[]- [Q. The 
latter equations can, in principle, be obtained from our 
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leading equation (Jllj) by neglecting statistical correla- 
tion, so as to average every factor independently. 

Using our effective-medium method, we have investi- 
gated the relaxation of charge carriers in band tails at 
T = 0. According to our results, energy relaxation is 
connected with dispersional transport. Even if we have 
no real diffusion in the system, since the temperature is 
zero, we have some spreading of the energy distribution 
with time. For a constant density of states the situation 
is completely equivalent to the one in the problem of the 
electron's motion in a disordered system, subjected to 
an electric field [ p^[ . The main difference between the 
above mentioned problem and the energy relaxation is in 
that, while in the former case the sites are distributed 
homogeneously in space, in the energy relaxation prob- 
lem the density of states is a decaying function of energy 
for most physical systems. In this situation we arrive 
at the picture that the particles of a package with lower 
energies move slower than particles with higher energy. 
This leads to the opposite tendency: at first, there is a 
slowing of the package spreading. The time dependence 
of the dispersion, being linear at the first stage of the 
evolution, is later slowing down. Later on, different pos- 
sibilities exist, dependend on the particular type of the 
energy dependence of the density of states. 

We performed a detailed investigation for the expo- 
nential density of states for two time regimes. If the 
variation of the velocity v in the energy space across 
the distribution is smaller than the velocity itself, that 
is if Lop- / A <C 1, we arrive at the situation of a pack- 
age, moving steadily down along the energy axis without 
any deformation. Dispersion becomes time independent 
(see Eq.([3l|)). However, since the particles are sinking 
down this condition becomes violated when time goes 
by. The parameter cup 2 / A, being small at the first stage 
of the evolution, is getting large. When this parameter is 
larger than one, the steady motion condition becomes vi- 
olated again. The package, being of gaussian form before, 
undergoes some restructuring to another, non-gaussian 
stable form, with its width lower than before by some 
numeric factor of the order of 1 (Eq.(90)). In our pic- 
tures, Figs it is clearly seen, that the package be- 
comes non-gaussian. This result remains valid until the 
very moment the quasi-elasticity conditions break down 
MA«1). 

For the exponential density of states we have found 
the packet to move as lnlni, i.e., its motion is strongly 
slowing down with time, the packet becomes, roughly 
speaking, almost stopped. This type of behaviour may 
be called "glassy", because of the overall time scale for 
the packet evolution (governed by the exponential func- 
tion of the large parameter A/cj), becomes huge. For 
example, when ui/A = 10 -4 , this time scale approaches 
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The main simplification used in our paper is the quasi- 
elastic approximation. This approximation relies on the 
smallness of the upper bound of the energy transfer to the 
phonon system by one hop. For localized electrons, this 
upper bound can be much smaller as the Debye energy 
of the host material, since not all phonons can interact 
with localized electrons equally well. For localized elec- 
trons the electron phonon coupling constant approaches 
zero, for phonons with wave-vector q > 2a. Thus high 
energetic phonons are less effective. Only phonons with 
energies u) < W£>2aa, where a is the lattice constant of 
the host material, are effective. Furthermore, in disor- 
dered systems the high energetic acustical phonons are 
localized, and thus need not contribute to transport by 
necessity. Nevertheless, the question whether the quasi- 
elastic approximation is applicable depends much on the 
material of interest. If, however, we compare our result 
for the time dependence of the mean energy (^2|) in an 
exponential density of states with those existing in the 
literature , we also find agreement for to/ A w 1 , which 
indicates that our results are, at least qualitatively, of 
wider validity. Unfortunately, to our knowledge, in the 
literature there are no further results on the width of the 
energy distribution available to compare with. 

From our point of view the main question remained 
open is how the results, obtained for the exponential den- 
sity of states, may be generalized for other types of energy 
dependencies. One can imagine, e.g., that for densities of 
states decaying with decreasing energy slower than expo- 
nential, the dispersion is growing, but under some sublin- 
ear law. Also, in the opposite case, in a density of states 
decaying faster than exponential, we would possibly have 
a dispersion tending to zero for large times. 

It should, however, also be mentioned that it is not 
completely clear under which situation zero temperature 
results can be applied to systems at finite temperature. 
In our approach we assumed that T = 0. However, in 
a real system, when the carriers are sinking down, the 
criterion for the temperature to be treated as zero, is 
violated at every finite value of T at some moment of 
time, since the contribution of hops to sites with higher 
energy values becomes more and more comparable with 
one of downward hops. Therefore, the consideration of 
temperature becomes vital. Work in this direction is in 
progress. 
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APPENDIX A: ON THE CONFIGURATION AVERAGE 



Here we derive the set of integral equations (|17|)-(|2lj). The same method has already been used in Ref. [ |19| . 
According to Eq.(|l0|), the configuration average of the electron density (n(p, s)) is given by: 



(n (p, a)) = J dpip Q (p x ) P ( P i,p) , (Al) 

where 

P(p',p)=r,(p')$(p',p) (A2) 
is the diffusion function. According to Eq. (p"l|) , the diffusion function is given by: 



p 



S n V S 



Diagrammatically this series can be represented as depicted in Figj^. In this picture every full dot represents a 
potential V, which depends on disorder via the structur factor 77. The expansion starts with a single structural factor 
77, depicted by the empty dot. Note that all disorder is comprised into the structure factor 77. 

In order to calculate the configuration average, we average the series ( |A3| ) term by term. The resulting expansion 
can be depicted in the usual form, as in impurity scattering problem (see Fig.||). One can see from the picture, 
that the set of all diagrams can be decomposed into two subsets, the set S containing all diagrams connected with 
the empty point and the set F containing all other diagrams. Owing to this decomposition configuration averaged 
diffusion function can be written as: 

(P(p', p)) = J d Pl S if/, Pl ) F ( Pl ,p) , (A4) 

from which Eq.(|l7j) follows. Note that F is nothing but ($). 

Again the class of diagrams contributing to F can be decomposed into reducible and irreducible diagrams. All 
irreducible diagrams can be comprised into a function II. Doing so, we obtain Eq. (p8|) . Note, that, although the 
introduction of the irreducible part II is parallel to the introduction of the self-energy it has to be stressed, that its 
physical interpretation is completely different. In the present context Tl(p' , p) is the true transition probability from p' 
to p. It is fl(p', p) that has to be calculated and compared with the experimental situation and not the bar transition 
probability, as suggested in many papers (see e.g. |3|). Diagrams, contributing to II are depicted in Fig ]3L The set 
of diagrams depicted in Fig. |^ can be generated by means of the propagator ITp, defined by the equation (|2C|), using 
Eq.©. 

The propagator lip can also be used to generate the irreducible block S. Diagrams contributing to S are depicted 
in Fig.|[ From the picture it follows, that S is given by Eq. (^l]) . 

APPENDIX B: TRANSPORT COEFFICIENTS 

According to Eq. (J35|), v(e, s) is given by: 

v(e,s) = \u?N (e)W (q = 0\e) = \u 2 N » ^cxp [- Pc ( £ , s )] / dR— =— L r ^ T , (Bl) 

2 2 J 1 + exp [2aR — p c (e, s)\ 

where: 

f(e,s)u = expp c (e,s) . (B2) 

We assume, that p c is large, so the integrals can be calculated in the limit p c — + 00. This assumption is justified 
a posteriory by explicit calculation of p c . In limit p c — > 00, the integrand reduces to a step function, so that the 
integrals can be calculated easily. Doing so, we find: 



Pc (g,s) 

2a 



vexp[-p c (e,s)] , (B3) 
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where S(d) is the solid angle and d is the spatial dimension. From the calculation and the result it is to be seen, that 
p c (e, s) is to be identified with the dimensionless critical hopping length at energy e. 

The same procedure is applied to the calculation of D{e,s), defined by Eq.(B4). Using the same approximations, 
we obtain: 



1 S{d) 

D{£ ' s) = 2d(dTY) 



oAf(e) 



Pc (E, S) 



2a 



d+2 



Thus D{e, s) and v{e, s) differ only in preexponential factors. 
We anticipate that for s = we obtain: 



Pc (e,0) 



2a 



2da 



vexp[-p c (e,s)] . 



1/d 



(B4) 



(B5) 



where a is a number and S(d) is the solid angle. A different form for the diffusion coefficient can be found in Ref. 
|pOf . The expression given in this paper differs from our expression essentially only in that oj is replaced by the tailing 
parameter. In looking on this difference one should, however, take into account, that the model used in that paper is 
different. Furthermore, in Ref. [^0| the transport coefficient are independent of s. It is therefore not completely clear 
to which range of time they apply. Clearly, they can not determine the evolution for all times. 



APPENDIX C: DERIVATION OF SELF-CONSISTENCY EQUATION 



Here we derive the self-consistency equation (f>7\). In order to work out explicitely the first order contribution to 
the self-consistency equation ( |56| ) we first need the first order correction to II. It is given by: 



n (1) (p',p;s) = J 



dpxdp-zdpiN {e 3 )w P3 (//, Pi) s) F (pi, p 2 ) w P3 (p 2 ,p;s) 



(CI) 



If we insert the expression for w P3 , we obtain: 



n (1) (pt,p;s) 



dpi 



[M (e) W {pi, p; s) [P (p, pi) - F {pi, W (pi , p; s) - JV ( ei ) W {p' , pi; s) [P { Pl , p) - P (p', p)] W (p, Pl ; s) } . (C2) 
Thus, at this stage the self-consistency equation takes the form: 

= J dqde/de! {ei - e) {W (e) W (0\sf, e; s) P {q\e, ei) W (q\e 1} e; s)-Af (e) W {q\e', s; s) P (q\e', ei) W {q\e u e; s) 
-N (ei) W (0|e', e 1 ;s)P {q\s u e) W {q\e, £l ;s)+M (ei) W {q\e> , e 1 ;s)P {q\e' , s) W {q\s, sr, a)} . (C3) 

In order to eliminate P, we have to replace P{q\e', e) = F{q\e' , e) — f{e, s)6{e' — e), where F is the effective medium 
approximation of the diffusion propagator. Let us first work out the local contribution. Replacing F by / * S, the first 
and the third term of Eq.(C3) cancel each other. The 4th term is zero, taken into account that the delta function of 
the effective medium is multiplied by its argument. Thus, when F is replaced by /(e, s)S{e' — e), only the 2nd term 
of Eq.(|C3|) survives. To simplify this term we take into account that 



W{ q \e; S ) = W{0\e, S )^{^^), 



2a 



(C4) 



where is a dimensionless function. Consequently, we obtain 
dqde' [e' - e) Af (e) / (e', s) \w {q\e', e; s) 

i! 



r/q 



W{q\e, S )\ —f{e,s)N{e) 



1 



2a 
Pc{e,s) 



M{e)f{e,s)W 2 {0\e,s) / d d x4> 2 {x) 



(C5) 
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Here terms proportional to uif 1 (e, s)df(e, s)/de <C 1 have been neglected. 

Let us now focus on the contribution of the regular part of the diffusio n p ropagator to the self-consistency equation. 



Owing to the step functions in F an W, the first and third term in Eq.(C3) are zero, when F is replaced by F. Thus 
we are left with: 

J dqde/d £l (et - e) j-W (e) W {q\et, e; s) F {q\et, £l ) W {q\e x , e; s) + j\f(e\_) W (q\e' , e l ;s)F {q\e', e) W {q\e, e l5 s) } . 

(C6) 



The range of integration in ( |Cq ) is determined by the step functions in F and W(q\s' , e; s) = 6{e' — e)9{ui — e' + 
e)W(q\e';s). Owing to these step functions, the energy integrations extend at most over intervals of length u>, so 
that the quasi-elastic approximation can again be applied to the effective-transition probabilities and the density of 
states. The diffusion propagator enetering Eq.( |C6| ), however, can not be dealt with in this way since for arbitrary q 
the derivatives of the diffusion propagator are not small as compared to the diffusion propagator itself. Therefore, 
another procedure is needed. To simplify this expression further we consider Eq.(^). If terms small with respect to 
ujN' /N and uiW'(q\e, s))/W(q\e, s) are neglected the function 

2a 

$(x\y',y)=u J 2 F(x— -\u>y' ,uy)W(ojy,s)Af(u>y), (C7) 

can be introduced, that satisfies the equation 

1 

8 $(x\y',y) = 6(y'-y)+ [ dyi[<Z>(x\y', y x + y)cf>(x) - $(x\y', y)0(O)]. (C8) 

o 



Then the expression ( JC6| ) can be cast into the form 

It/ 1 l—j/i 

wW(uy)( , 2a , ) d [ d d x<t> 2 {x)[- [ dy' f d yi y'$(x\y' + -, Vl + -) + [ [ dy'y'$(x\y' + -, -)] (C9) 
Pc(uy,s) J J J u u J J uj uj 

oo oo 



For s — Eq. flCq ) does not contain any physical parameter. It only leads to the determination of the function 

®o(x\y',y) = Q(y' - y)®o(x\y' - y), 

which satisfies the equation 

l 

o = s( y '-y) + J d yi [<t>o(x\y',yi + y)<t>(x)-®o(My',y)m}- (cio) 

o 

Provided, we restrict our consideration to small frequencies in deriving the self-consistency equation we only need to 
take into account the linear contribution of the function $ with respect to s/(W(ui y )J\f(ujy)uj). Then, using again the 
smallness of the variation of the effective transition probabilities and the density of states with respect to changes of 
energy over intervals of length uj, <!> can be approximated as 

y 



$(x\y', y) = $ (x\y' - y) - — — — — / d yi <S> (x\y - yi)^ (x\ yi - y). (Cll) 

W {u)y)N {ujy)uo 



y 



Using this expression the self-consistency equation takes the form (pjl) , where the coefficients a and b are given by 



l 



Jd d x<j> 2 (x)f(l-yf$ (x\y) 
1 o 



a ~2 fd d x^(x) ' (C12) 

_ 1 J d d x^(x) Jp 1 dy{\ - yf \l dy&Mv ~ yM^Vi) ™ 
2 JdV(x) ■ 1 ' 
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APPENDIX D: CALCULATION OF THE ENERGY RELAXATION SPEED FOR CONSTANT DOS 



To calculate the integral (|72j), we again use Eq. (|6q) to change the integration variable from s to y — v/vq. Further- 
more, for convinience, we introduce the parameter r = fit. Doing so, we obtain: 



dE (t) 
dt 



-v F(t) , 



where F(t) is given by: 



F ( T )=[ p- l + iny exp (rylny) 
J c 2m \ny 



After an integration by parts we obtain: 

dF 



1 



dr 



dy 



1 



— = / — exp (ry lny) = S (r) 



r J c 2m' 



(Dl) 



(D2) 



(D3) 



To calculate F(t), Eq. (D3) should be supplemented with the condition F(+oo) = 1. 

To evaluate S(t), it is convenient to use the integration contour Im(ylny) = 0, or, introducing polar coordinates 
y = r exp (i9): 



r = exp (—6 cot 9) 



Substituting Eq.flDj) into Eq.flpJ), we have: 

** dS 



S(t) = 



-e cote 



(D4) 



(D5) 



The value of the integral may be estimated by the saddle-point method, looking for the maxima of the expression 
\n.V (8) — tV (9). The stationary point equation is: 



V (9) 



1 



V{9) 



0. 



As r> 1, the stationary point is 8 S =0, and the asymptote of Eq.(|D5|) is: 

S(t) ps (27reT) _1/2 exp(-r/e) . 



(D6) 



(D7) 



On the other hand, if r <C 1, we have two stationary points ±0 S , V {±9 S ) = 1/t, 9 s — tt — 5, and for S <C 1 the 
following equation may be obtained: 

7T /7T \ 1 

the solution for which at small r is: ir/5 = W(e/r) w In (e/r) — lnln(e/r) .The asymptotics of S (r) turns out to be: 



S(t) 



pn(e/r)]" 



(D8) 



APPENDIX E: THE SADDLE-POINT APPROXIMATION 



For energy dependent density of states it turns out to be difficult to obtain explicite expressions for the time 
dependence of the energy distribution function. Here, according to Eq.(f4(i|), the time-dependence has to be calculated 
from the equation: 



+ 200 



P(s,s ,t) 



ds 

2niv (e, s) 



exp 



st 



de' 
v (e', s) 



(El) 
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where v(s, s) is given by Eq. ( |66| ) or ([38]). Let us try to integrate over s, using the saddle-point approximation. The 
saddle-point position So(e, £o,t) may be obtained from the equation: 



ds' 



v(e')g(e',s ) 



so 



de' dg (e, so) 

v (e') g 2 (e' , s ) <9s 



(E2) 



where <? = s)/v(e, 0) = exp W(s/f2(e)) was introduced. The expression for the diffusion propagator in the saddle- 
point approximation may be written as: 



P L (e,eo,t) 



■ exp 



ds' 



dg (e, so) 



°J v (e')g 2 (e\s ) ds 



(E3) 



where: 



D(s,s ,s ) 



d 2 f de'n(e') In g{e',s ) 



ds 2 



2v (e') 



> 0. 



(E4) 



Assuming so/fi <C 1, we have: g (e, so) ~ 1 + so/Cl — s 2 l /2Q 2 

HQ 



so 



t-T(e,e Q ) 
2D (e,e ) 

and, finally, the diffusion propagator: 

P L (e,£o,*) 



D(e,e ) 



V (£') (£') 



eo 

f de 

-, T (e, Eq) — / — - 

) J v {e') 



(E5) 



^4ttL> (e,e ) v ( £ ) 



exp 



(t-T(g,e )r 
4D (e,e ) 



(E6) 



At a given time t the exponent in the above distribution function is maximal at e = e m , where e m is given by the 
condition t — T(e m ,£o)i or: 



de' 
vo (£') ' 

Expanding the expression (Eq) around s = e m , we have gaussian distribution: 



(E7) 



Pl (e,eo,t) 



V (£m,Eo) V27T 



exp 



(g - £ m (£p,t))' 
2cr 2 (£ m ,£ ) 



(E8) 



with the dispersion: 



o 1 (e™, £o) = 2wq (e m ) £> (e m , £o) = 2w 2 (e m ) 



cfe' 



(E9) 



Note, that according to Eq.(E7), the dispersion is time dependent. If the density of states Af(e), and, consequently 
0(e) and vo(e), are strongly varying functions of e, the integrals can be further simplified. In this case we have: 



de' 
He 1 ) 



dv (e ) 


-l 


dv (e m ) 


de 




dp 



(E10) 



ds' 



V (£') (£') 



dfr (£q)»(£q)] 

de 



~ d\v (e m )n(e m )] ~ 



-1 



(Ell) 
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de' 



v (e') n (e') 



dv (eg) Q (go) 
de 



efo (e m ) (e m ) 



(E12) 



The applicability condition for the saddle-point method is: 

SI (e') In g(e',s ) 



D d (e,e ,s ) » 



de' 



2v (e') 



(E13) 



Note, that this condition restricts the applicablity of the saddle-point approximation to times, that are not too large. 
Also, from Eqs.(E5) one can see, that this approximation is invalid when e is close enough to the initial point £o, — 
here |so| becomes large, and it is not possible to expand over Sq. So, the initial stage of the evolution, when the entire 
distribution is concentrated near £q, has to be investigated separately too. 



APPENDIX F: CALCULATION OF THE MOMENTA FOR LARGE TIMES 



Here we present details on the calculation of the momenta of the distribution function for large times and for the 
exponential density of states. In our calculation we take into account that the density of states decays with decreasing 
energy. 

Let us change the integration variable in (|88|), p — > z = W (spe p ), p = W (s _1 ze z ). As p 3> 1, we have, introducing 
p s = W(s- 1 ), s = p7 x exp(-p s ): 

p = W (p s ze ps+z ) » p s + z + lnz , 
as \z\ <C |yO s | and |lnz| <^ \p s \. Both inequalities holds within actual region of integration. We can write now: 



Xul(s) 



-- I dz 



lnp s 



z + In z 



Ps 



d 

a~ ex P 

oz 



1 f, 1 
1 + -z 



(Fl) 



Something more about the approximation we have chosen: every extra power of z/p s means an extra power of \Q, 
which can be neglected. An extra power of z means the multiplier ~ \f^p s . All our considerations are valid if uip s <C 1 
only. Therefore, we have the reason to omit extra powers of z/p s ~ \[Q in the following, but to keep some extra 
powers of z~ x ~ l/y/c5p s 3> yfiL. 
Thus we have, for example: 



Xi (s) = — In p s H / dz (z + In z) — exp 

s sp s J dz 

o 

ln(Vu;p s ) lny 



= — Inp 

s 



dy \y 



l 

upl 
d 



z[l + -z 



y/ujp s Vtups I dy 



exp 



T_ 

2 



VtOPs 



(F2) 



where 



Vujp s 



Within three leading orders on the parameter \j\fuip s we have: 



Xil(s) 



In p s 



ix'jj 1 In ( v ~'l>- 
2 s sp s 



ln2- 7 \ 1 



sp s 



which can be written, using p a + lnp s = ln(l/s) also as 
Xil(s) 



7T 1 

8o) sp 2 s 



_lnln J — - + J — —,b = V2e ^ 2 , 

s s V 2 s y 8Ll) spj. 



(F3) 



(F4) 
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where 7 is the Euler's constant. In the time representation we have: 



2-kis 



I TTOJ I TT f 

V ~2 + V &5 Jc 



ds 



c 27risln 2 & 



At large t, by/ujt 3> 1, the first integral, after the substitution: st = y, may be expanded as: 



L 



ds 
c 2iris 
1 



In by/ujt Jc Zniy 



dy y, 
e" In y 



dy 
2-Kiy 



In 



In b\fujt — In y e v = In In b\fujt 
dy 



2 In by/tit Jc 2, *iy 



3 y In 2 



y- ... 



In In b\fuit + ■ 



7 



+ 



^ 2 /6- 7 2 



Inb^Cdf 21n 2 feV^i 

As for the second integral, one may restrict oneself with the leading term only. As a result, we have: 

T 1 



Xi (t) — — In In b\[u)t - 
Let us consider now the second momentum 



+ 



8w In 2 by/tit 



+ ...,b= V2e- 1+ ^ 2 . 



OO 

X2l{s) = J ^ln 2 p^ L (p;.s)^(ln 2 p) . 



(F5) 



(F6) 



Again, we change the integration variable to z = W (spe p ), write p 3> 1 as p = p s + z + In z, and set the lower limit of 
the integration to be zero. Restricting ourselves with the lowest order in z/p s ~ we set lnp = \np s + {z + \nz)/ p Sl 
and: 



dp , dz d 

—<i>L = — ^- cx p 

p s oz 



1 

up 2 s 



z\l + \z 



Taking into account: 



we can write: 



Xil(s) 



X2l{s) = s X \ L {s) 



--In p s 

s 



z + In 2 



+ In z 



z + In 2 



(F7) 



Then we have: 



z + hiz 



/ dz (z + lnz) 2 -^- exp 



spt 
s 



dz 



upt 



z[l + -z 



dy [y 



In (\/u)p s ) my \ y 2 y 



\fuip s V^Ps 



exp 



2 \fuips 



Expanding this expressions in power series on l/y/ujp s , and keeping three leading powers, it may be easily obtained: 
'z + lnz\ 2 \_2(I> \/27rcj In (v^Ps) In 2 + 7 y/2nui In 2 (Vtip s ) 



+ 



SPs 



( _ ln2_7\ In (V^p s ) 

V 2 ) sA 



2 sp s + sp 2 s 
2" 



In 2 + 7 + 

24 ' 



In 2 - 7 



1 



,2 ' 
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^From Eq. (F3) we have: 



: + hi : \ / Trio 1 \n(^/u}p s ) / ln2 — 7\ 1 
2 s sp s \ 2 J sp s 



So, from Eq. (|F7|) we obtain: 



X2l(s) = xIl(s)+ 2-- +(l-ln2) 



However, we want to calculate central momentum: 

=X2(*)-Xl(*): 

or, in the Laplace representation: 

M2l(s) =X2i(s)- (Xi) L 0), 

+ioo+0 



\J~2-KLO I IT" 7T 

24 + 2 ~ 



7T 1 

8w sp 2 



^) l ( S )ee jdte st xl(t)= J gl m ( Sl ) m(s - Sl ) 



-ioc+Q 



Taking the expression for xiz, from Eq. (F4), we have 



" V+0 dsi lnln^lnln^| 6 ^ ™ 

^ 7 f^ + ~ lnln^— + — 

Zm s\ (s — s\) s s Is 

zoo+0 



+ioo+0 



ds\ 



In In 



b\fu> 



2ni 



si(s-si)ln 2 -^ 2s In 



2 by/qj 



Let us consider the first integral. After changing the integration variable, s\ — sx, we have: 

1, 2l b^ 2, bV* f dx^fi-Mi-*)/^^)) 

- m m 1 — m ' 



dx 



s s s 
In ( 1 - In 



2?ri 



x/ In (^)) In (l - ln(l - x)/ In (^)) 



2ttz 



x(l — x) 



(F8) 



(F9) 



(F10) 



(Fll) 



Closing the integration contour to the left, one can easily prove, that the first intehral in the above expression is zero. 
The second one at small s may be expanded into powers of 1/ Ins. In the leading order we have: 



b\fu) 



— In 2 In 

s s 



1 



dx lnxln(l — x) 



j n 2 fey 7 " J c 2-iri x(l — x) 



Performing integration by part in the integral, and then deforming the integration contour to make it run along the 
lower and upper shores of the cut (-co, 0) of the integration function: 



dx lnxln(l — x) 1 

~2 



c 2ni x(l — x) 



dx 2 
- — In x 

q 2m dx 1 — x 



d ln(l-x) 1 f , ln(l + a;) - 1 _ . 2/ 

-— / dx^rz — Imln (x + iO) 



2tt 



= dx 



[ln(l 



1] In: 



(1 + xf 



The evaluation of the second integral in Eq. (Fll) gives: 

-2 — V + O 



1 



s In 



2 b\fuj 



s In 



4 b^O 



24 



As a result, one can write: 



(xl) L (s) = s X 2 1L ( S ) + 



6sp 2 s 



Thus, one have the second central momentum in the Laplace representation to be: 



M2L » = 2 " o - + 1 " In 2 Y- + - - — - 1 )-= 

2 s sp s V 2 8 / spj. 



and the time dependence of the second central momentum: 



Ma(*) = (2-|)w + (l + In2) 



V27 



hit 



(F12) 



(F13) 



(F14) 



Expression for the third central momentum p,^ (t) (^lj) may be obtained quite analogously. 



APPENDIX G: ASYMPTOTIC FORM OF THE DISTRIBUTION FUNCTION AT LARGE TIMES 



One can write the large times distribution function ( p7|) , after the replacement of the integration variable s — > z 
W (s/?e p )as follows: 



(M) = J^-(M) 



(Gl) 



-\-ioc-\-5 



i>(e,t)=peP 



dz d 
— 2( 1 + z )— ex P 



2ttiz 



de 



-z + t-e z ~ p - — z \ l + -z 
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bjp* 



(G2) 



— ioo-\-S 

One can try to calculate the latter integral within the saddle point approximation. The saddle point equation is: 



f(z) = -l + t 



z+l , 



-p 



- —s{z+l) = 0. 

Uip 



Assuming \z\ <C wp 2 (remind that top 2 ^> 1), we can neglect the third term in the above equation, and to obtain: 



z c + 1 = W 



As we shall see later, the actual values of parameters p, t, corresponding to the body of the distribution function, obeys 
the inequality pe p /t ^> 1, and, therefore, z c ^> 1. But this leads to the saddle point parameter, f"'{z c )/ [/"(z c )] 3 ^ 2 = 

\(z c + 3)/(z c + 2) 3 / 2 | w 1, to be of the order of 1, instead of being much less. Let us note, however, that the 
function in the exponent in Eq. (GS) may be represented as a sum of two terms, of which the first one, — z + 
tz/pexp [z — p), ensures the function in the integral to have a maximum at z = z c of the width ~ 1, while the second 

one, (tip 2 ) 1 z (1 + z/2), may be supposed to be almost constant within this peak, — it varies essentially at scales 
Az ~ V&p S> 1 within the actual interval of system's parameters. Thus, one can replace that second term in the 
exponent by its value at some zq = z c + a, \a\ ~ 1. After this, the integration may be easily performed, returning 
back to the initial integration variable s: 



t/,(M) 



d 

-Ptt exp 
dp 



ojp* 



;W(s peP) 



1 + -W (s pe») 



ds 
litis 1 



— 



where the last integral is equal to simply t, and, taking into account Zq S> 1, Sq = Zo/P ex P(^o — P) ~ c Ai where 
c ~ 1 will be obtained later. Performing differentiation with respect to t, see Eq. (G2), and neglecting zo/ojp 2 <C 1, 
we obtain finally: 
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-4^-^<^)[ i+ MH} 



(M) 



To obtain c, let us calculate Xi(i) = e(t) with the distribution function ( |G3| ). We have: 



(G3) 



= ~ J -^4>(p,t)fop = J dz lnp(z) ^exp 



W I -ze' 



where the replacement of the integration variable p — > z = W {cpe p /t) was done. Assuming /)> 1, one can write, using 
W (a;) « a; — In x as |x| ^> 1, that p(z) « W (t/c) + z + \nz = p t + z + lnz (the condition z ~ "C Pt = W (i/c) was 

used). In the leading order on \/ui, one can replace p{z) within the exponent with p t , and In p{z) — * lnp t + (z + ln z)j p t . 
The limits of integration by z are to be set (0, +oo). Then we have: 



DC 

Xi{t) = -lnp t - J dz- 



+ In z d 
Pt dz 



exp 



1 ( 1 

ujpf V 2 



(G4) 



This latter integral, analogously to one in Eq. ([F]), may be evaluated as an expansion on small parameter \j^fwp: 

XX (*) = - (V2e-^^f) - Jf +. • • (G5) 



Comparing Eqs.(|G5|) and 189), one can see, that c = exp (—7). 
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